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Abstract 

We consider a long, semiflexible polymer, with persistence length P and contour length L, 
fluctuating in a narrow cylindrical channel of diameter D. In the regime D <C P <C L the 
free energy of confinement AF and the length of the channel Ru occupied by the polymer are 
given by Odijk's relations AF/R\\ = AokBTP^^D- 2 / 3 and R\\ = L [l - a Q (D/P) 2 / 3 ] , where A Q 
and a are dimensionless amplitudes. Using a simulation algorithm inspired by PERM (Pruned 
Enriched Rosenbluth Method), which yields results for very long polymers, we determine A Q and 
a and the analogous amplitudes for a channel with a rectangular cross section. For a semiflexible 
polymer confined to the surface of a cylinder, the corresponding amplitudes are derived with an 
exact analytic approach. The results are relevant for interpreting experiments on biopolymers in 
microchannels or microfluidic devices. 
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I. INTRODUCTION 



Microfluidic devices provide new possibilities for studying biological polymers such as 
DNA, actin filaments, and microtubules. Since the persistence lengths of biological polymers 
are typically tens of nanometers or larger, their behavior in confinement, as in nano- or 
microchannels, is different from that of flexible synthetic macromolecules. 

In this paper we consider the equilibrium statistics of a semiflexible polymer or worm-like 
chain with persistence length P and contour length L in a channel of diameter D. In the 
regime D <^ P <^ L, corresponding to a long, tightly confined polymer, Odijk [lj showed 
that the free energy of confinement AF, i.e. the work required to reversibly insert the 
polymer in the channel, and the length of the channel R\\ occupied by the polymer are given 
by 
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For a channel with a rectangular cross section with edges D x and D y , 



^-.M J_ + J_) ( 3 ) 



/ n 2 / 3 _i_ n 2 / 3 

R\\ = L\l-a a p V f v ) . (4) 



Here A Q , a Q , A n , and an are dimensionless universal numbers, which do not depend on P, 
D, D Xl and D y . 

Making use of advances in the manipulation of single polymers, recent experiments have 
begun to approach the Odijk regime D <^ P <^ L. In the experiments of Reisner et al. [2] 
on single DNA molecules with persistence length P of about 50 nm and contours lengths 
L of around 20 fim or larger, the condition P <^ L is well satisfied, and the dimensions 
D x = 30 nm, D y = 40 nm of the narrowest channels are moderately smaller than P. In the 
experiments of Koster et al. [3] on actin filaments with persistence length of about 20 fi in 
microchannels with diameters down to 1 or 2 /im, D <^ P for the narrowest channels, and 
the longest contour lengths L considered of around 50 fim are about 2 to 3 times P. For 
an experiment in which DNA is confined by a grooved substrate instead of a channel, see 



Hochrein et al. 



a. 



For interpreting such experiments it is important to know the numerical values of the 
dimensionless amplitudes in Eqs. (G])-(IIJ). Solving an integral equation numerically which 
arises in an exact analytic approach, Burkhardt [5] found 

A n = 1.1036 , (5) 

and from simulations Bicout and Burkhardt [6] obtained 

A u = 1.108 ± 0.013 , A Q = 2.375 ± 0.013 . (6) 

Other estimates from simulations, compatible with these values but with larger error bars, 
are given in Refs. jl, 8, 9], and related results for a helical polymer in a cylindrical channel 
inRef. |loj|. 

Although the free energy amplitudes A u and A Q are known with good precision, compa- 
rable estimates of the extension amplitudes a n , a have not been available. Thus, we have 
determined a a and a Q from simulations and obtained new, more precise estimates of A u 
and A QJ as described in this paper. 

The surface of a channel can be prepared so that biopolymers are adsorbed. For ex- 
ample, naturally anionic DNA strands are adsorbed on a surf 
membranes and have a high lateral mobility on the surface 
action between the surface and the biopolymer lowers the free energy barrier for insertion 
of a macromolecule in a narrow channel. With this as motivation we also consider the free 
energy and extension of a semiflexible polymer confined to the surface of a cylinder. 



ace coated with cationic lipid 
. The attractive inter- 



The theoretical framework for our calculations is outlined in Section II. In Section 
simulation algorithm inspired by PERM (Pruned Enriched Rosenbluth Method) jl5 . 



16 



I a 



is described, which enables us to consider polymers two or more orders of magnitude longer 
than in the simulations of Ref. [6|, on which the results (JHJ) are based. Our estimates of 
the amplitudes An, A Q , a a , and a Q are given in Section III. In Section IV we consider a 
semiflexible polmer confined to the surface of a cylinder with diameter D. In the regime 
D <^ P <^ L the confinement free energy and extension are also given by Eqs. (CQ) and 
{2D, but with different amplitudes As, as- These amplitudes are calculated with an exact 
analytic approach. In the concluding Section V we compare the results of Sections III and 
IV with predictions for a polymer confined by an effective parabolic potential. 
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II. THEORETICAL FRAMEWORK 



In the worm- like chain model of a semiflexible polymer, the bending energy is given by 



Here f is the unit vector tangent to the polymer contour, s is the arc length, and k is the 
bending rigidity, related to the persistence length by P = K/kgT. In the regime D <^ P 
L, backfolding of the polymer and excluded volume effects are negligible. In typical polymer 
configurations the tangent vector is nearly parallel to the symmetry axis of the channel. 
The configurations correspond to single valued functions r(t), where (x,y,t) are Cartesian 
coordinates (see Fig. 1), and r — (x,y) specifies the transverse displacement of the polymer 
from the symmetry axis or t axis of the channel. Since \v\ <C 1, where v — dr/dt 1 the 
bending energy ([?]) simplifies to 

„, « r L , (d 2 f\ 2 

and the length of the channel R\\ occupied by the polymer and the contour length L are 
related by 

L= [ * dt(l + v 2 ) 1 ' 2 ^R\\ + \ [ ll dtv 2 . (9) 
Jo 2 J 

In accordance with Eq. (JSJ), the partition function of a polymer with position and slope 
r , at t = and r, v at t is given by the path integral 



Z(r, V] r , v ; t) = J D 2 r exp ~ dt 



(10) 



where r is restricted to the interior of the channel. It satisfies the Fokker-Planck type 
differential equation 

d 1 



ft + * ' Vr " 2P V ' ) ^ ^ = ° ' 

with the initial condition Z(r, v; r , iT ; 0) = S(r — r )5(v — v Q ). 

The boundary condition at a "hard" channel wall follows from the fact that that dis- 
continuities in the slope of the polymer cost an infinite bending energy and are suppressed. 
Thus, as r approaches the channel wall, Z(r J v;r ,Vo]t) vanishes for n • v > 0, but not for 
n • v < 0, where n is normal to the wall and directed toward the interior of the channel [5]. 
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Our reason for denoting the Cartesian coordinates by (x,y,t) instead of (x,y,z) is ex- 
plained in Fig. 1. Each polymer configuration r{t) may be interpreted as the position of a 
randomly accelerated particle in two dimensions, plotted as a function of the time t. The 
polymer partition function ( II (J corresponds to the propagator or probability density for 
propagation from initial position and velocity ro,^o to f,v in a time t. From the Boltz- 
mann factor in Eq. (fTUI) one sees that the acceleration of the particle at each instant is an 
independent, Gaussian-distributed random variable, with 

$~ = #), (#)> = o, (vMvAt')) = S f s(t - 1') . (12) 

Since the polymer partition function vanishes at a hard wall for n • v > 0, the propagator 
for the randomly accelerated particle vanishes if the particle is reflected toward the interior 
of the two dimensional domain representing the channel cross section. Thus, the hard wall 
in the polymer problem corresponds to an absorbing boundary for the randomly accelerated 
particle. 

For large t the partion function (fTUI) decays as 

Z{r, v- f , *T ; t) « Mr, *#o(r , -v^e'^^ , oo , (13) 

where Eq(P,D) is the smallest eigenvalue of the t independent Fokker-Planck equation 

1 



-E(P,D) + v-V r -2pVl 



i[)(r,v) = 0. (14) 



Together with the definition 



6XP| k B T)~Z{oo) U5j 



of the free energy of confinement, Eq. (1130 implies 



^(s.l) 



k B TR l{ W fl )-W»> = #^ (16) 

for a channel with a circular cross section, in agreement with Odijk's relation (OQ). Here we 
have used the scaling relation Eq{P 1 D) = (2P)~ 1 / 3 D~ 2 / 3 E'q (|, 1), which is readily derived by 
rewriting Eq. (HM in terms of the dimensionless variables r' = D _1 r, t' = (2P)~ 1 ^ 3 D~ 2 ^ 3 t J 
v' = (2P) 1 /3 jD -V3^ 

For a channel with a rectangular cross section with edges D X1 D y , the partition function 
in Eq. (TTU1) has the product form Z(r, v; r , v ; t) — Z(x, v x ; x , v x0 ; t)Z(y, v y ; y 0l v y0 ; t). This 



5 



is the origin of the sum of independent x and y contributions in Eqs. (|3| and (J1J). The 
solutions to Eq. ( TH1) also have the separable form Tp(r,v) = i/j(x,v x ,)ilj(y,Vy) J implying 
E°(P, D x , D y ) = E l (P, D x ) + E l (P, D y ). Here e[(P, D x ) is the smallest eigenvalue of the 
equation 

r f) 1 f) 2 1 

il>{x,v x ) = 0. (17) 



, did' 



dx 2P dv 2 

on the one-dimensional interval — \D X < x < \D X , with boundary condition ip(—^D x , v x ) = 
ijj(\D x ,-v x ) = for v x > 0. The scaling relation E l (P,D x ) = (2P)- 1 / 3 D~ 2/3 E l (^l), is 
readily derived by rewriting Eq. (TT71) in terms of the dimensionless variables x' — D~ l x 1 
t' = (2P)- 1 / 3 D x 2/ \ v' x = (2Py/ 3 D x 1/3 v x . 

From Eqs. (OP), Q, (TTBI) . and the results of the preceding paragraph, we obtain 

A = 2- 1 / 3 ^(i,l), ^ D =2- 1 / 3 4(i,l). (18) 

To obtain comparable expressions for the amplitudes a Qj a n , we begin by comparing Eqs. 
(ED and (|9]), which imply 

«o = 2 J (v%b • (19) 
Here (v 2 )° PD is the average value of v 2 along an infinitely long, tightly confined polymer in 
a channel with a circular cross section. In terms of the ground state eigenfunction ipo(r,v) 
of Eq. (HID with eigenvalue E°(P, D), 

^2y = I d2r I d2y g 2 jW 1 v)Mr, ~v) ^ 
j d 2 r j d 2 v Tjjo(r,v)iljQ(r, —v) 

Expressing Eq. (1191) and its analog for the rectangular cross section in terms of the dimen- 
sionless variables introduced below Eqs. ( TT6l) and ( fTTl) . we obtain 

a„ = 2- 5 / 3 (i;' 2 >l 1 , a u = 2-^ {v V)\ , (21) 

2' J ' 2' 1 

Equations (I18j) and (1211) play a central role in our work, allowing us to determine the free 
energy and extension amplitudes from simulations with P = ~ and D — D x — D y — 1. 

III. SIMULATIONS 
A. Algorithm 



To determine A u and an from simulations, we generate a large number N of configura- 
tions x'(t') of a polymer with persistence length P = = in the unbounded two-dimensional 
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space (x',t'). Here x' and t' are the dimensionless coordinates introduced below Eq. (JT 
The configurations are generated with same Boltzmann weight as in Eq. OlOp . but in two 
rather than three spatial dimensions. All of the configurations have the same initial position 
and slope x' = v' = at t' = 0. Each configuration is "grown" until it leaves the interval 
— ~ < x' < | for the first time. From this information we calculate the fraction Q(t') of the 
N configurations which have not yet left the interval at t'. 

From Eq. ( TT5|) we see that Q(t') = exp(— AF/ksT), where AF is the free energy of 
confinement of a polymer with one end fixed, as described above, which extends a distance 
i! down a two-dimensional channel of width 1. According to Eq. (1 13ft and the discussion 
below Eq. (fTTl) . Q(t') decays as 

Q(t') ~ e-^id' 1 )*' (22) 

for large t'. To estimate A a , we fit the Q(t') extracted from the simulations with the 
exponential form ( 1221) for large t' to obtain £q(|, 1) and then use Eq. (JTSJ) . 
The polymer configurations are generated with the algorithm 

/A 3 \ 1/2 / \ 

X n+1 = X 'n + < A n+l + ( ~ ^ J \S n +l + ^3 T n+1 ) , (23) 

v' n+1 = v' n + (2A n+1 ) 1/2 r n+1 , (24) 

introduced in Ref. and also used in Refs. [s, 12]. Here x' n is the position of the 

polymer at point t' n , and A n+1 = t' n+1 — t' n is the length step. The quantities r n and s n are 
independent, Gaussian random numbers with (r n ) = (s n ) = and (r^) = (s^) = 1. 



As discussed in Refs. 
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this algorithm generates polymer configurations consis- 
tent with the Boltzmann weight (1 10 J) in free space, i.e., in the absence of boundaries. An 
advantage of the algorithm is that in free space there is no length- step error. The length 
step A n+ i need not be small. For good efficiency we use a fairly large step when x' n is well 
inside the interval — \<x'<\. Near x' = ±| a smaller step is needed in order to accurately 
determine the value of t' at which the configuration leaves the interval for the first time, and 
hence Q(t'). As in Ref. [6], we choose 

A n+1 = 10" 1 (| - |x|) + 10" 5 , (25) 

which varies from 0.05 at x' = to 10~ 5 at x' — ±|. Further reduction of the length step 
had no significant effect on our estimates. 
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In estimating (v 12 ) J in Eq. (12TT) to determine one should only use the subset of the 

2 ' 

configurations, generated as described above, which lie entirely within the interval — § < 
x' < |, i.e., within the channel. For each of the configurations C in the subset we average 
v l2 1 over the entire length of each of the configurations in the subset, using 

n 

Here t' = A n+1 is the total length of the channel occupied by configuration C, and 
(v' 2 ) n+ i is the equilibrium value of v' 2 for a semiflexible polymer with endpoints (x' n , v' n ) and 
{x' n+1 ,v' n+1 ) at t' n = and t' n+1: respectively, averaged over all intermediate t' . This quantity 
is readily calculated from the free space partition function or propagator and is given by 

(v' 2 )n+l = TE^n+l + | (x'n+l ~ X n) + ]| ( V n+1 ~ \ V 'n+l V 'n + V n) 

~\ K+l - <) (<+l + V'n) A nll ■ (27) 

Having calculated (v l2 )c for each configuration in the subset this way, we average the results 

over all the configurations in the subset to obtain an estimate of (v' 2 ) J and, using Eq. ( |2T1) . 

2 ,i 

the corresponding value of a D . 



B. Enrichment Procedure 



The quantity {v' 2 )\ in Eq. (I2T1 ) is the average value of v 2 for a semiflexible polymer 
of infinite length in a channel. We found it necessary to go to lengths t' of around 100 
to estimate (v' 2 } J v free of finite- length effects, to 3 significant figures. However, it is not 



,/2\ 

2 ' 

feasible to generate configurations this long, which lie entirely within the channel, without 



modifying the steps outlined in the preceding paragraphs. From Eqs. © and ( |2T1) . i?o(f , 1) 

is close to 1.390. Thus, according to Eq. (I22l) . the probability that a configuration of length 

t'=100, generated as described above, never leaves the channel — | < x' < |, is about 
e -139 w 1Q -6i 

To generate a large, statistically useful number of configurations lying entirely in the 
channel, we used an enrichment procedure inspired by PERM (Pruned Enriched Rosenbluth 

ations of a wide variety of 



Method) \m 



Il7j , which has been successfully applied in simu 
systems, including flexible, self-avoiding polymers in channels jl8. 
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We begin by generating a large number Nq of configurations as described above. Let N\ 
be the number of these configurations which have not yet left the channel at t' = r. We 
make n copies of each of these configurations and then, with the algorithm of the preceding 
Subsection, continue each of the nNi configurations past t' = r. Let N 2 be the number 
of these configurations which have not yet left the channel at t' = 2r. Again we make n 
copies and then continue the nN 2 configurations past if — 2r. At if — 3r, 4r, . . . the same 
procedure is followed. 

To estimate A a using Eqs. (1T81) and (1221) . we need to calculate the probability Q(t'), 
defined above Eq. (1221) . that a configuration, generated as in the preceding Subsection, has 
not yet left the channel at if . To obtain this probability, it is useful to think of copying all 

the preceding configurations, including the number of initial configurations, at t' = r, 2r 

Thus, 

Ni. 

Q(0) = 1, Q{kr) = - j ^r forfc = l,2,... (28) 

Our results for a semiflexible polymer in two dimensions were obtained with ]Vq = 1.8 x 
10 7 , r = 1, and n = 4. These values of n and r were chosen so that N k slowly decreases 
with increasing k. From Eqs. @, (flgjl . (1221) . and (I25|> . one finds N k ~ n k e~ E ^^ kT NQ = 
(0.996) fc A^o • 

We have also calculated the number of families N[ am to which the N k configurations that 
remain in the channel up to if = kr belong. Two configurations are said to belong to the 
same family if they coincide in the interval < if < r, i.e., if their most remote ancestor is 
the same. By definition Nl am = N u but for larger k, N{ &m < N k , since several of the A/fc 
configurations may belong to the same family. According to our simulation data A^| am also 
decays as n k e -E (~,i)kr _ (o.996) fc . For sufficiently large k all N k configurations belong to a 
single family. 

To determine a D , we evaluate (v'^}[ v as outlined above in the paragraph containing 

2' 1 

Eqs. ( f26l) and (1271) . for those N k configurations which remain in the channel up to if = kr, 
estimate the limiting value for large t', and then use Eq. (12T1) . 

The simulations of a polymer in a channel with a circular cross section of diameter D 
are very similar. In terms of the dimensionless Cartesian coordinates (V, y f , t') introduced 
below Eq. ( TT6l) . the channel has radius ~. Beginning with x' = y' Q = v' x0 = v' y0 = t' Q = 0, we 
generate the sequence [af n ,y' n ,if n ) with the algorithm (I251) - (I2^I) and corresponding equations 
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with x replaced by y. In analogy with Eq. (1251) the length step is 

A n+1 = 10" 1 



|-« 2 + ^ 2 ) 1/2 ]+ 10 " 5 - ( 29 ) 
Each configuration is grown until it leaves the circular domain [x a + y' 2 ) 1 ^ 2 < \. Again 
we begin with N configurations and at t' — r, 2r, . . . make n copies of the A r 1 , N 2 , . . . , 
configurations which have not yet left the circular domain. Our results were obtained with 
N = 4.7 x 10 7 , r = 1.009, and n — 20. As in the two-dimensional case these parameters 
were chosen so that Nj~ and A| am decay rather slowly, as n k e - E o{\^) kT — (0.997) fc , where we 
have used Eqs. JED, Q, (EE}, and our result for A Q in Eq. (1301) . 

To estimate A Q , we calculate Q(t') for integer t' using Eq. (1251) . fit the results with the 
exponential form (1221) . but with £? °(|, 1) in place of Eq(^,1), and then use Eq. (fT8l) . To 
estimate a , we evaluate (v /2 )^ , as described in connection with Eqs. ( |26l) and ( 1771 ). for 

2' 

those Nk configurations which remain in the channel up to t' = kr and then use Eq. ( 12 lj) . 
C. Results 

In Fig. 2, lnQ(t'), as determined from Eq. (l28i) . is shown for t' — 0, r, 2r, . . . for a 
polymer on a two-dimensional strip (upper curve) and in a three-dimensional channel with 
a circular cross section (lower curve). According to Eqs. (1221) and ( TTB1) the upper and lower 
curves have slope Eq(^, 1) = 2 1 / 3 A D and i?o(|, 1) = 2 1 ^ 3 A QJ respectively, for large i! . From 
the best fit to the slope, we obtain 

A u = 1.1038 ± 0.0006 , A = 2.3565 ± 0.0004 . (30) 

The uncertainty was estimated by determining A u and A Q in each of the intervals 100 < 
t' < 200, 200 <t'< 300, ... , 900 < t' < 1000 and quoting a value somewhat larger than the 
width of the corresponding distribution. 

The estimate for A n is in extremely good agreement with the result in Eq. (JSj), obtained 
jy solving an integral equation numerically that determines A a in an exact analytic approach 
5]. The new estimates for A u and A Q in Eq. (I30|) have smaller error bars than the earlier 



estimates [6] shown in Eq. (|6j), which are based on simulations of much shorter polymer 
chains. 

In Fig. 3 the t' dependence of (v' 2 ) J and (v /2 )f is shown. The averages are based 

2' 1 ' 2' 1 

on the configurations which remain in the channel from the starting point up to i! . For if 
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greater than roughly 100 the curves are consistent, within statistical fluctuations, with the 
constant values (v' 2 } J 1 = 0.2901 ±0.0003 and (v ,2 )\ 1 = 0.5400±0.0004. Substituting these 
values in Eq. (12TT) . we obtain 

an = 0.09137 ± 0.00007 , a Q = 0.1701 ± 0.0001 . (31) 

For a polymer on a two-dimensional strip, the probability distribution P((v /2 )) of the 
quantity (v' 2 ) is shown for representative values of t' in Fig. 4. The distribution was 
determined from the results for {v' 2 )c 1 where the index C labels the configurations that 
remain in the channel up to length t', and (v l2 )c is the average value of v' 2 along configuration 
C from the starting point up to t', calculated as in Eqs. ( |26l) and (1271) . The distributions in 
Fig. 4, are approximately Gaussian, and the half width or standard deviation w 1 shown in 
Fig. 5, decreases in good agreement with the t /_1 / 2 law expected for statistically independent 
contributions. Results similar to those in Figs. 4 and 5 were also obtained for a polymer in 
a channel with a circular cross section. 

Since (v' 2 )c determines the contour length L' of configuration C via Eq. (J9|, the curves in 
Fig. 4 may be interpreted as distributions of the contour length L' for fixed i! . Presumably 
the distribution of t' for fixed V ', i.e., the end-to-end distribution for a polymer of fixed 
contour length [20], is very similar. 

IV. SEMIFLEXIBLE POLYMER CONFINED TO A CYLINDRICAL SURFACE 

In this Section we consider a semiflexible polymer constrained to lie on a cylindrical 
surface with a circular cross section. As mentioned in the introduction, this is an obvious 
model for a semiflexible polymer adsorbed on a channel wall. We analyze the case in which 
only configurations that leave the polymer in contact with the cylindrical surface are allowed 
and the equilibrium statistics is determined by the Boltzmann weight e"^ fcBT , where TL is 
the bending energy (J7J) . The bending energy is clearly minimized if the polymer configuration 
is a straight line parallel to the symmetry axis of the channel. In the limit D <C P <C L the 
free energy of confinement and the extension again are given by Eqs. (CQ) and (J2j), as shown 
below, but with different amplitudes As, czs-, which we determine. 

In the regime D <C P L 1 typical configurations of a semiflexible polymer correspond 
to single valued functions r(t) = r{t)R 1 where r specifies the transverse displacement of 
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the polymer from the symmetry axis or t axis of the channel. Here (r, 9 1 t) are cylindrical 
coordinates, and f , i are the corresponding unit vectors. Keeping in mind that dr/d0 — 0, 
dO/dO = — r, and that \RdO/dt\ <C 1, one finds that the bending energy {7J) takes the form 



H ~2 



dt 



d 2 e 



dt 2 



1 



\ [R — 



R 2 



dd 



dt 



(32) 



As discussed below Eq. (fTTI) . the polymer configuration r(t) may be interpreted as the 
position of a randomly accelerated particle in the (x, y) plane, plotted as a function of time. 
A polymer confined to the surface of a cylinder corresponds to a particle moving on a circle 
of radius R. The first and second terms in the integrand in Eq. ( |32l) are the squares of the 
tangential and centripetal accelerations, respectively. 

Equation (1321) and the definitions x = RO — \D0 and v — \D dO/dt lead to the partition 
function 



Z(x - x , v, v ,i) = J Dxexp i ~ J dt (^^j 



dx 

^D 2 [~di 



(33) 



and the Fokker-Plack equation 

d d 2P 



1 & 



h V 1 V — 

dt dx D 2 2P dv 2 



Z(x — Xq, V, Vq, t) — 



(34) 



Disregarding the position of the polymer endpoint, we integrate Eq. ( |3ll) over x from — oo 
to oo. This yields the Schrodinger equation 



d 2P 
di + D- 2V 



1 



Z(v,v ,t) = 



(35) 



2Pdv 2 

Equation ( |35l) also follows directly from the path integral Z(v,Vo,t) = 
j Dv exp j — |P J* dt [(dv/dt) 2 + (4/D 2 )v A ] |, which has the same Boltzmann factor as 
in Eq. (l33l) . but expressed in terms of v rather than x. 

Beginning with Eq. (1351) and following the steps that led from Eq. (fTTI) to Eqs. (fT8l) and 
( I2TT ). we obtain 

A £ = 2-^E s {\, 1) , a s = 2- 5/3 <«' 2 >f t , (36) 

2 ' 

Here Eq(~, 1) is the ground state energy of the Schrodinger equation with a quartic potential 

d 2 



-E° + v' 4 - 



dv' 2 



i/>(v') = 



(37) 
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written in terms of the dimensionless variables introduced below Eq. ( TT7I) . ipo(v') = ipo(-v') 
is the wave function of the ground state, and {v /2 )f is the quantum mechanical expectation 

2 > L 

value of v' 2 in the ground state. Equation ( |36l) allows us to calculate the free energy and 
extension amplitudes from simulations with P = ~ and D = 1. 



The ground state energy of Eq. (1371 ). determined numerically by Bender et al. 2l|] and 
Voros |22j . is given by Ejf (|, 1) = 1.060 362 09. Solving the Schrodinger equation (|37|) 
numerically for this energy using Mathematical we obtain (v /2 )f = 0.362 023. Substitution 

2 ' 

of these values in Eq. (1361) yields the amplitudes 



As = 0.84161 , a s = 0.11403 



(38) 



V. CONCLUDING REMARKS 



Using a PERM inspired simulation algorithm, we have determined the amplitudes A n , 
A , and a Q for a semiflexible polymer in a channel with an estimated error of less 
than a tenth of a percent. We hope the results will be useful in analyzing experiments. 
For a polymer confined to the surface of a cylinder, we have calculated the corresponding 
amplitudes As, &s exactly to five significant figures with an analytical approach. These 
latter results may be used as a benchmark for testing simulation algorithms. 

A common approximation in studies of semiflexible polymers in channels is to replace the 
hard wall interaction by an effective parabolic potential 2Q , l23| . In concluding, we use this 
approximation to relate the free energy and extension amplitudes A and a and compare the 
relations with our simulation results. 

The partition function of a polymer which is tightly confined along the t axis by a 
parabolic potential energy per unit length u — \ ksTbr 2 is given by the path integral 



Z(f } v;r 0} v ;t) = / D 2 r exp 



1 



dt 



It can be calculated exactly |20|, |23| and yields 



\dt 2 



+ br 



/46 

f(b, P) - k B T 
(v 2 ) = (46P 3 )" 1/4 



1/4 



(39) 

(40) 
(41) 



where / = —hsT lim^oo t 1 In Z is the free energy per unit length, and (v 2 ) is evaluated in 
the same limit t — > oo. 
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For a parabolic potential the confinement free energy per unit length is Af = f(b, P) — 
f(0,P) = f(b,P). To obtain a prediction for a , we choose b so that this Af exactly 
reproduces expression (pQ) for the free energy of confinement in a channel with a circular 
cross section, evaluate (v 2 ) in Eq. for this 6, and then substitute the result in Eq. (TTU|) . 
This and a similar calculation for a channel with a square cross section yield 

a n A n = - , a A = - . (42) 
o z 

n 

Ubbink and Odijk [24] argue that the parabolic confining potential is an artifice, and 
that the confinement free energy (OD) of the polymer in a channel should not be identified 



with the full free energy / in Eq 
/ - bdf/db = §/. This leads to [25 




but with the configurational part /config = f ~ ( u ) 



3 3 

oluAu = — , a A Q = - , (43) 

which differs from Eq. (1421) by an extra factor of | on the right hand side. Our simulation 
results in Eqs. (|30l). (13T1) . and ( 1551) yield the products 

a n A n = 0.1009 ± 0.0002 , a A = 0.4008 ± 0.0003 , a s A s = 0.095969 , (44) 



which lie in between the predictions ( 1421) and (1431) but in every case closer to (1431) . Since the 
surface of a cylinder is two-dimensional, the product as As should be compared with the 
predictions for a u A u in Eqs. and (1151) . which apply in both two and three dimensions. 
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FIG. 1: The curve may be interpreted as a tightly confined semiflexible polymer in a channel with a 
circular cross section or as the trajectory r(t), plotted as a function of t, of a randomly accelerated 
particle moving in two dimensions, which has not yet left a circular domain. 
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FIG. 2: lnQ(t') vs. t' for a polymer on a two-dimensional strip (upper curve) and for a polymer 
in a three-dimensional channel with a circular cross section (lower curve). 
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FIG. 3: {v'^)\ vs. t' for a polymer on a two-dimensional strip (lower curve) and (v vs. t' 

2' 1 2, 

for a polymer in a three-dimensional channel with a circular cross section (upper curve). 
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FIG. 4: Distribution P((v' 2 }) for a polymer on a two-dimensional strip, as defined below Eq. (j31 
The curves correspond, from bottom to top, to t' — 100, 225, 400, 625, and 900. 
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FIG. 5: Double logarithmic plot of the half width or standard deviation w of the distribution 
P((v' 2 )) (see Fig. 4) as a function of t' . The straight line corresponds to w = kt'~ l f 2 , where k is a 
constant. 
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